Time dependent couplings and crossover length scales in non-equilibrium surface roughening 
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We show that time dependent couplings may lead to nontrivial scaling properties of the surface fluctuations 
of the asymptotic regime in non-equilibrium kinetic roughening models . Three typical situations are studied. 
In the case of a crossover between two different rough regimes, the time-dependent coupling may result in 
anomalous scaling for scales above the crossover length. In a different setting, for a crossover from a rough to 
either a flat or damping regime, the time dependent crossover length may conspire to produce a rough surface, 
despite the most relevant term tends to flatten the surface. In addition, our analysis sheds light into an existing 
debate in the problem of spontaneous imbibition, where time dependent couplings naturally arise in theoretical 
models and experiments. 

PACS numbers: 81.15.Aa,05.40.-a,64.60.Ht,05.70.Ln 



in 



i 

C 

O 

o 



> 

in 

m 

^6 
o 
i> 
o 



X 



The theory of surface growth has applications to phenom- 
ena that cover a wide range of length scales from nanometers 
to millimeters, including the growth of thin films from an in- 
coming flux of atoms fllL Si , fluid imbibition in porous me- 
dia |4j], and propagation of fracture cracks (5[], among many 
others. It is often observed that surfaces kinetically roughen 
and become scale-invariant. 

In kinetic roughening details of the interactions are largely 
irrelevant in the mathematical description of the critical prop- 
erties of the surface at long wavelengths, akin to critical point 
phenomena. Therefore, in the hydrodynamic limit, one can 
describe surface growth in d + 1 dimensions by the stochastic 
equation 



•r?(x,t), 



(1) 



where h(x, t) is the height of the interface at substrate position 
x and time t. The functional Q(Vh) depends on the specific 
model and should satisfy all the symmetries and conservation 
laws. The external noise ry(x, t) describes the random driving 
forces acting on the surface, for instance the influx of particles 
in a deposition processes. 

We shall be discussing here surface growth models with 
local coupling among degrees of freedom, i.e, growth equa- 
tions that include only terms that depend on derivatives of 
the height, T = V/i. In the spirit of Ginzburg-Landau- 
Wilson theory of critical phenomena the functional Q(T) is 
constructed as the leading-order expansion in powers of the 
argument T, its derivatives, and combinations thereof. One 
explicitly avoids to include terms that are incompatible with 
the symmetries of the problem. The corresponding expansion 
takes the form of a sum, Q(T) = J^i Wi<Pi(x, t), where a;, are 
the coupling constants and <£>i (x, t) are the local operators 
in the terminology of the renormalization group (RG). Local 
operators correspond to the combinations of degrees of free- 
dom and its derivatives, which typically include surface dif- 
fusion V 2 /i, curvature diffusion V 4 /i, and in general higher- 
order diffusion terms V 2m h. Also nonlinear terms usually ap- 
pear, (V/i) 2n , V 2 ™(V/i) 2m , and so on. The asymptotic long 



wavelength limit is governed by the most relevant terms in 
the RG sense. Higher-order terms are nonetheless important 
in describing crossover effects before the truly asymptotic be- 
havior is reached. This leads to the existence of crossover 
length scales and characteristic times at which one can ob- 
serve the true asymptotic scaling behavior of the system. In 
the case of time independent coupling parameters the analysis 
of crossover length scales is relatively simple and very well 
studied . However, little is known about less common 

situations in which the couplings depend explicitly on time. 
Examples include the problem of a stable phase growing at 
the expense of a metastable phase |g], and spontaneous imbi- 
bition Jfjj], among others. Indeed, the problem of spontaneous 
imbibition has been a subject of great interest in the last few 
years ABBS [IS [ill]- Different theoretical approaches have 
arrived at the conclusion that, for small deviations around of 
the mean position H(t), the fluid-fluid interface (in Fourier 
space) is given by 



dh k (t) 
dt 



= -aK\k\k 2 h k - H{t)\k\h k + K\k\ff k , (2) 



where a and K are the surface tension and permeability con- 
stants, respectively, and the average position follows Wash- 
burn's law H(t) = y/ H 2 + 2at and arises from mass con- 
servation. There is an interesting debate regarding the role of 
this time-dependent coupling in the scalin g ob served in both, 
simulations |7l [KJ ll 111 and experiments 1U2L 1 1 3 O. This 
problem has largely motivated our study of the scaling of the 
surface fluctuations in systems where couplings depend ex- 
plicitly on time. 

In this Letter we study the interplay between a crossover 
length scale growing in time and the dynamic correlation 
length characterizing the kinetic roughening process. We find 
nontrivial scaling properties, including anomalous roughen- 
ing il5lll6ll . The long time limit scaling behavior strongly de- 
pends on the nature of the phases that the dynamic crossover 
separates. We focus on three typical cases that cover the most 
important situations one can find: (i) Crossover between two 
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different rough regimes, (ii) crossover from a rough to a flat 
regime, and (iii) the existence a damping term, where scale- 
invariant fluctuations are damped over a certain length scale 
that varies with time. We exemplify our general scaling anal- 
ysis with simple model systems in which dimensional analysis 
gives the exact exponents. In order to gain analytical under- 
standing we restrict ourselves to linear model examples that 
allow exact computation of the critical exponents. Results are 
compared with numerical simulations. 

Crossover between two different rough regimes.- Con- 
sider a growing surface described by Eq. ([T). Invariance un- 
der translation along the growth and substrate directions as 
well as invariance in the election of the time origin rule out 
an explicit dependence of Q on h, x and t. These symmetry 
requirements alone lead to scale invariant growth in a generic 
fashion |l|] . Let us consider the two most relevant terms in the 
leading-order expansion of Q{ Vh) in the hydrodynamic limit, 
so that we have Q = oji&i(x, t) + uiu&u(x, t) + higher-order 
terms. Let the $1 term be more relevant (in the RG sense) 
than $n- We first discuss the case in which Eq. (Q~|i exhibits 
crossover between two different rough regimes. This means 
ujn$ii(x, t) is relevant at short scales, while uji&i(x,t) be- 
comes the most dominant in the long wavelength limit. Hence 
one expects to find an early times (short scales) regime with 
scaling exponents oS^ and z^ iv> that crosses over to the true 
asymptotic regime, governed by the operator $1, with expo- 
nents aW and Dimensional analysis indicates that there 
exits a crossover length £ x ~ (wn/wi) 1 '', for some exponent 
q. This is the typical length above which the most relevant 
term wi$i takes over. This is likely the most common situa- 
tion of crossover behavior in surface growth. 

The question we wish to address here is how this picture is 
modified when the most relevant operator's coupling, u>\, de- 
creases in time. The most interesting situation naturally arises 
for a power-law decay of the coupling u>\{t) ~ i -7 , so that 
the term u>i(t)$>i gets effectively less relevant for longer times 
as compared with the oju^ii term. In this case the crossover 
length diverges in time as a power-law, £ x (t) ~ i 7 / g , and 
the interplay with the other relevant scale in the problem - 
namely, correlation length £(t) ~ i 1 / 2 , will give rise to an 
interesting behavior. On scales smaller than £ x (t) the less rel- 
evant term II dominates, while we expect a crossover to the 
asymptotic regime I for ^> t x (t). The caveat is now that 
the observable scaling regimes crucially depend on the expo- 
nent of the coupling 7. If 7 < q/z^ the crossover does take 
place as described above. On the contrary, for 7 > q/z^ the 
crossover length grows at an exponentially faster rate than the 
correlation length. Therefore, in this case the time-dependent 
coupling effectively takes the crossover length to scales much 
larger than those that can be correlated by the dynamics at 
any finite time. As a consequence the system will never cross 
over to regime I, and only the "less" relevant operator II gov- 
erns the scaling regime observable for arbitrarily large system 
sizes. 

We now study a simple 1 + 1 dimensional model that exem- 
plifies the crossover behavior between two different scaling 



regimes above discussed. We consider ([T]) with $1 = V 2 /i 
and $11 = ~V 4 h, describing diffusive coupling. Then we 
have the growth equation 

dh 

= v (t)V 2 h-V 4 h + r){x,t), (3) 

at 

where we have rescaled all the couplings but v(t) to unity. The 
noise is delta-correlated, (r)(x, i)rj(x' ,t')) = 25(x — x')S(t — 
t'). Higher-order diffusion terms, V 2n /i, might be present, 
but are irrelevant for the scaling behavior. Actually, even the 
V 4 /i term is irrelevant as compared with V 2 h if the couplings 
are independent of time. Consider now a decaying coupling 
v(t) ~ t~ 1 [17]. Balancing the two gradient terms in Eq. (01 
one finds that they become comparable at the typical crossover 
scale £ x (t) ~ t 1 ! 2 , with the exponent q = 2 in our previous 
analysis. 

Dimensional analysis tells us that the kinetic roughening 
process is governed by the fourth derivative term (that we 
label as II) on scales smaller than £ x , while the diffusion 
term should dominate on much larger scales (that we label as 
regime I). Below the crossover, for <C i x (t), Eq. (01 has 
a correlation length £(t) ~ i 1 / 2 ' , where z' 11 ) = 4, given by 
the V 4 /i term dynamics. Crossover to the asymptotic regime I 
takes place when the correlation length reaches the crossover 
length t x l z ~ Vl 2 . Following our previous analysis we 
expect that the regime I can only be reached if 7 < 2/z^ u \ 
i.e, 7 < 1/2. On the contrary, for 7 > 1/2 the crossover is 
wiped out and the less relevant operator II (in RG sense) does 
govern the scaling in the long wavelengths limit as well. This 
analysis is in excellent agreement with simulations as shown 
in Figure 1 . 

Let us now study the scaling properties in the case 7 < 1/2, 
where one expects two different regimes with different crit- 
ical exponents. Smaller scales, regime II, are simply gov- 
erned by the dynamics of the curvature diffusion term V 4 /i 
and we find the well-known results z^ 11 ) = 4, = 3/2, 
and the local roughness exponent a?^ = 1 fl5l [liill . Above 
the crossover, £(t) ^> £ x , the system is expected to be in 
regime I and a scaling analysis gives the exact global ex- 
ponents a® = (1 + 7)/[2(l - 7)], = 2/(1 - 7), 
= (1 + 7)/4, and the scaling relation a (I) = z (I) /3 (I) 
is fulfilled. Figure Q] shows that this scaling analysis is in ex- 
cellent agreement with a numerical integration of the model 
©. 

However, even for 7 < 1/2 the asymptotic scaling be- 
havior is nontrivial, since the time-dependent coupling v{t) 
leads to anomalous scaling of the local surface fluctuations. 
This is more clearly proved in Fourier space in terms of 
the structure factor (or spectral power spectrum), S(k, t) — 
(h(k,t)h(—k,t)), where h(k,t) is the Fourier transform 
of the surface in a system of lateral size L, h(k,t) = 
L^ 1 / 2 ^2 x [h(k, t) — h(t)] exp(ikx). ^,From the growth equa- 
tion and neglecting the V 4 /i term we obtain 

S(k,t) - fc" 2 t 7 s(fc 2 i 1 ™ 7 ), (4) 
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FIG. 1 : (Color online) Crossover from rough to rough. Numerical 
integration of Eq. ® using a time step At — 10~ 2 and c = 10~ 3 
in a system size L — 1024. All results were averaged over 200 
realizations, (a) Global width of the interface, computed for dif- 
ferent 7 values. The theoretical prediction for the growth exponent 
P m = (1+7) /4 for 7 < 1/2 fully a grees with simulations. For any 
7 > 1/2 we obtain /3 = 3/8, as expected. The curves are vertically 
shifted for clarity, (b) Power spectrum for 7 = 0.3 evaluated at times 
10, 10 2 , 10 3 , and 10 4 . Anomalous scaling with a local roughness ex- 
ponent aifo = 0.5 is observed. The inset shows a data collapse of 
the local width data for a (I) = 0.93(1) and z (I) = 2.86(4) to be 
compared with the predicted values = 0.928 and = 2.857 
for 7 = 0.3. 



where the scaling function s(u) has the asymptotes, s(u) ~ u 
for u -C 1, and s(u) ~ constant for u ^> 1 [18]. The 
temporal shift in the power spectrum, f, implies anoma- 
lous scaling of the local fluctuations. The theory of anoma- 
lous scaling lfl5l [lrjll tells us that the scaling relation 7 = 
2(a (I) - a-ioJ/z® must be fulfilled. We then find that, ac- 
cording to the spectral density the scaling in the asymp- 
totic regime (regime I) is anomalous in this case with a local 
roughness exponent = 1/2. Note that the local rough- 
ness exponent does not depend on the coupling exponent 7 
and self-affinity, a® = a[P c = 1/2, is recovered whenever 
the coupling is time independent (7 = 0). A comparison with 
numerical simulations for 7 = 0.3 is shown in FigureQ] 

Crossover from a rough to aflat regime.- Another inter- 
esting situation occurs when the most relevant operator hap- 
pens to lead to a flat surface. In this case the crossover 
is expected to take place from an early times (short scales) 



rough regime II, with roughness exponent cr 11 ) > to an 
asymptotic flat regime, where cr 1 ' = 0. Again, a time de- 
pendent crossover length leads to nontrivial scaling proper- 
ties in this case. Below the crossover we expect a rough 
phase dominated by the <!>n term with exponents cr 11 ) > 
and z^ Il \ Above the crossover, however, surface correlations 
cannot evolve any longer since the operator $1 does not am- 
plify surface fluctuations (note that a™ = and z^ = 0). 
This leads to an asymptotic roughness exponent identical to 
that in the early regime. Moreover, it is worth to stress here 
that the dynamical length scale ~ t x l z 1 becomes a constant 
above the crossover (note that zW — 0). Therefore, above 
the crossover the only relevant scale left in the problem scales 
as ~ Vl q . This must be identified with the correlation length 
£(t) above the crossover. We therefore conclude that the time- 
dependent coupling leads to an asymptotic regime with expo- 
nents a — and z = g/7, instead of those that would 
naively arise from the operator I, a® = 2® = 0. It is re- 
markable that, despite the most relevant term tends to flat- 
ten the surface, the time-dependent coupling conspires to pro- 
duce a rough surface. The roughness exponent is inherited 
from the early time (short scale) phase. Meanwhile, the time- 
dependent crossover length becomes the only relevant scale 
above the crossover, which is to be associated with the dy- 
namical exponent z of the system in the long times regime. 

To exemplify this general scaling analysis we now study 
a simple 1 + 1 dimensional model system exhibiting a 
crossover between rough regime on short scales to a flat sur- 
face regime in the long wavelengths limit. Let us consider 
same growth model as in (O, but the noise is now conserved, 
(r} c (x, t)r] c {^, t')) = -2V 2 5(x - x')8{t - t'). Note that the 
surface diffusion operator V 2 h leads to flat surface fluctua- 
tions, a® = = 0, in the presence of conserved noise. 

On small scales the surface is expected to be dominated by 
the less relevant V 4 /i term, so we obtain z^ = 4 and a^- 11 ' — 
a[^} = 1/2 in the early times (short scales) regime. Above 
the crossover the V 2 h term governs the surface fluctuations 
andwefindzW = 2/(1-7) and or^ = ~(l-37)/[2(l-7)] 
by simple power-counting. Therefore, for values < 7 < 1/3 
we have the desired situation of a crossover from a rough to 
a flat regime. According to our scaling argument we expect 
scaling in the asymptotic regime with a roughness exponent 
a = = 1/2 and a dynamical exponent z = 2/7, in- 

stead of the trivial values a' 1 ' = and = that would 
naively correspond to the V 2 /i operator with conserved noise. 
A comparison with a numerical integration of the model for a 
time-dependent coupling with exponents < 7 < 1/3 shows 
an excellent agreement with this analysis (see Figure 2). 

Crossover to a damped regime.- Finally, we briefly dis- 
cuss the case in which surface fluctuations become damped 
over a certain scale. This means that, for small deviations 
around the mean surface height, the terms —h, —h 2 , ... are 
to be included in the growth equation. These terms break 
the h — > h + c symmetry and therefore also break scale- 
invariance. Let us consider a growth model with the leading- 
order expansion Q = —uj(t)h(x,t) + oju^h(x, t)+ higher- 
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FIG. 2: (Color online) Crossover from rough to flat. Numerical inte- 
gration of Eq. (0 with conserved noise in a system size L — 1024. 
Results were averaged over 200 realizations, (a) Global width for 
different 7. The growth exponent was found to fit with the theoreti- 
cal value P = 7/4. The curves are vertically shifted for clarity, (b) 
Power spectrum for 7 = 0.3 at times 10 3 , 10 4 , 10 5 , and 5 x 10 s . It 
shows a roughness exponent a = 0.48(3). The inset shows a data 
collapse of the local width data using a — 0.51(2) and z — 6.63(4), 
which agrees with the scaling prediction 2 = 2/7 f° r 7 — 0.3. 



order terms. In such a way that the operator $n is less relevant 
than the —h damping term, but it would lead to scale-invariant 
behavior in the absence of damping. Following a scaling anal- 
ysis similar to that of the previous sections we find that for 
z^ < q/j the crossover to the damping regime will never 
take place. So, in all respects, the damping term is irrele- 
vant and the surface will exhibit the scaling behavior corre- 
sponding to the operator II. On the contrary, if > q/j the 
crossover occurs when £(t) ~ £ x (t). The asymptotic regime 
in this case does exhibit scaling but with the nontrivial expo- 
nents a — a' 11 ) and z = q/j, for the same reasons as in the 
previously discussed rough-to-flat case. Simulation results of 
the simple model (not shown), dfh = —u>(t)h + V 2 h + 77, are 
also in excellent agreement with this scaling analysis. 

Conclusions.- We have shown that the presence of time- 
dependent coupling in non-equilibrium surface roughening 
has important implications in the scaling properties of the 
asymptotic regime. We have focused on local models that in- 
clude the two most relevant terms separated by a crossover 
length scale growing in time, £ x (t) ~ Vl q , which conspires 
with the correlation length of the system £(t) ~ t x / z to give 
highly nontrivial scaling properties. In the case of a crossover 
between two rough regimes, the surface may exhibit anoma- 
lous roughening, directly related to the value of 7 ^ 0. On 



the other hand, in the case of a crossover from a rough to ei- 
ther a flat or damping regime, the dynamical crossover length 
£ x (t) ~ t 1 ^ is the on ly relevant scale above the crossover 
and this immediately leads to an asymptotic rough regime 
with a dynamic exponent z — q/7. A rough regime appears 
despite the most relevant term tends to flatten the surface. Re- 
markably, this is precisely the numerical result found in re- 
cent studies in spontaneous imbibition 0, OH EH, where 
the interface growth is described by Eq. (0. Our results show 
that these numerical results can be understood in the wider 
context of kinetic roughening in systems with time-dependent 
couplings. 
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